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Abstract 



^ . 

qh. We represent the method of dynamic diffusion for the approximate solution 

of Shredinger equation with decoherence. Decorence shows as the divergency 
of exact solution from the dynamics of diffusion swarm, which arises when 
. the total number of real particles grows. The method of dnamic diffusion 

cannot be reduced to the solution of differential equations, in contrast to 
• Bohm's quantum hydrodynamics, hence the direct computer simulation is the 

. mandatory step of its development. We compare the dynamic diffusion with 

l/^ I the explicit finite differences scheme for Shredinger equation for some standard 

' one dimension cases. 

1 Introduction 

, Experiments on quantum computers (see |T0]) show the necessity to develop quan- 

tum theory in the area of many particles for the investigation of decoherence prob- 
lem. It is closely connected with the general theory of complex systems and pro- 
cesses, in which quantum effects play the key role. For a systemn with n real particles 
we can write Hamiltonian of interaction in the form of general formula, not in the 
form of a certain matrix, because such a matrix has the dimensionality exp{n), e.g., 
for the accessible values n it is inaccessible due its huge size. The possibility to write 
the explicit formulas fro Hamiltinian H factually, gives us also the way to investi- 
gate fast quantum algorithms without launching them in reality. For example, when 
solving the search problem, e.g., finding the solution xtar of equation f{x) = 1 for 
the Boolean function f of n variables, we know neither /, nor Xtar, hovewer we can 

convince that f -v/y time iteration of Grover operator G gives us the state \xtar), 
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provided we start from the state with the equal amphtudes (see Here is the 
force of cartesian mathematics, e.g., of the formulas. 

Hovewer, if we apply this result practically, e.g., to the solution of the certain 
equation f{x) = 1, we immediately meet the principal difficulty that is the prac- 
tical realization of the operator G, which acts in the space of = 2" dimensions. 
Therefore, to find the root Xtar of the equation we must deal physically with Hilbert 
space of exponential dimensionality. At the same time for the establishing of the 
fact of convergency of the algorithm it is not necessary to do it and only formulas 
refiecting general properties of G, as its unitarity, and its concrete form G = I^Ixtar 
are sufficient. 

Experiments on quantum computers shows the existence of serious limitations to 
its scalability, which are typically called decoherence. Decoherence should be treated 
as the divergency of the real evolution of quantum system from the exact solution of 
its Shredinger equation. On the other hand, in the standard apparatus of quantum 
theory there are no intermediate forms of evolution between unitary dynamics and 
measurements and the last need the external subject, which has rights to measure 
(or to observe) a state of the small system at hand. We thus treat the observer not 
only as somethig (somebody) more complex than our system, but as the subject, 
which we cannot include to the investigated system in principal. 

The standard formalism of quantum theory thus establishes the border between 
"simple" objects, to which we can apply quantum mechanics and "complex objects", 
to which quantum theory is non applicable. This conclusion is new because it follows 
from our non ability to build a scalable quantum computer at least at this moment. 
Fundamental role of decoherence forces us to look anew to the attempts of L. de 
Brogile and D. Bohm to build the explicit urn scheme for quantum probability. More 
exactl we can name it the building of the " constructive version" of quantum theory, 
in which it can be generalized to "complex objects" (see [B|). 

The method of collective behavior rests on the idea of D. Bohm about the rep- 
resentation of a quantum particle as an ensemble of classical point wise particles, 
which we call its samples and the ensemble - the swarm. Each sample has not only its 
coordinates, but its own velocity. The genesis of this idea gets rize to L. de Brogile, 
who acquired the fundamental sense to the wave corresponding to the real particle, 
which has the dynamical characteristic - the velocity. The impulse of the swarm 
in a given point becomes then connected with the phase S/h of the wave function 
by the relation VS* = p. Resting on this correspondence E. Shredinger derived his 
famous equation. However, the swarm representation of quantum particle meets the 
following fundamental difficulty. In Shredinger equation occurs the wave function 
\E', connected with the density of samples by the equality p = and we thus have 
no direct connection between the dynamics and p. To make sense the application 
of the swarm we have to require that its density will be proportional to the density 
of theh probability of the occurrence of the particle in a given point. Then samples 
must be presenred in the time, e.g. each of sample must have its own history. The 
dynamics of the swarm consisting of behaviors of all samples must approximate the 
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exact solution of Shredinger equation during sufficiently large time frame. 



2 Bohm approach 

The problem of conservation of samples individuality has been solved by D. Bohm, 
who proposed the swarm of classical particles, which are driven by the external 
potential V = V{f) and the so called quantum potential Q = — where \& = 

R exp{iS/h). Shredinger equation with the external potential V is equal to two real 
equations of the form: 

l + v-O) =0, 

where the ffist is Hamilton- Jacoby equation for the characteristic function S of 
a sample, which is defined by the equations VS = p, dS/dt = E, where p = 
p{f) is its impulse, p(f) = B?{f) - the density of swarm, and the second equation 
shows the property of continuosity of a swarm. This approach is called quantum 
hydrodynamics, because a quantum particle is here represented by a swarm of its 
classical samples, which obeys the law of classical dynamics but their potential 
U = V + Q has the quantum summand Q, which depends on the density of the 
swarm in this point in the time and space. 

The algorithm modelling quantum dynamics by hydrodynamical approach at 
each step does the following. 

• Fulfills the free flight of all samples. 

• Change the speed of each sample accordingly to the rule p(t + At) = p{t) —VU. 

We see that the computation of the density of swarm at each step is very impor- 
tant, because it is required for the computation of ?7, due to the addition in the form 
of quantum pseudo potential Q. We can consider that samples merely determin the 
nodes of division of configuration space, which is the base of the finite differences 
method for the solution of ([1]); the value of quantum hydrodynamics is thus limited 
by the possible acceleration of the solution of Shredinger equation. How valuable 
is this adventange of quantum hydrodynamics? We consider the case of many real 
particles, for example a few tens of them. The density of swarm will be then the 
density of swarm in i?^" space, where n is the total number of particles. To conserve 
the initial accuracy of the approximation to 

pw^i^wr (2) 

for the solution ^ of Shredinger equation, we need to keep constant linear spatial 
step, which results in the exponential grows of the total number of cells in the 
configuration space. To be able to compute derivatives from the density (we need 
three sequential) we must guarantee that each cell is full of samples that erects 
the same difficulty as the direct solution of Shredinger equation for n particles. The 
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single profit from the quantum hydrodynamics consists in the possibihty to optimize 
the grid for the finite differences scheme; this adventage is interesting but it does 
not resolve our main problem. 

3 Swarm with dynamic diffusion 

Explicit consideration of mechanism for the acceleration of samples requires the 
satisfaction of conservation laws. If we are going to refuse from the knowledge of 
density in the vicinity of the point at hand, we must ensure these laws for separate 
samples, which means the refusal from such things as the quantum pseudo potential. 
First of all we define this mechanism and write the equation on the string of samples 
through the border, which this mechanism induces; then we show, that Shredinger 
equation can be reduced to this equation if we fix the grain of the spatial resolution. 
It turns that our mechanism of dynamic diffusion critically depends on the choise of 
this grain that causes some peculiarities of our approach. However, these peculiari- 
ties are inavoidable price for adventages we obtain with dynamic diffusion for many 
particles. 

3.1 Mechanism of bond creation and annihilation 

We will define the evolution of swarm, which we call the dynamic diffusion and which 
will consist of steps, each of duration At. But here the simple scheme from the 
previous section is not sufficient, because the swarm will substantially transformed, 
hence we define the step of evolution anew. 

At each step we will have not samples, but symplexes. A symplex S of an order 
j (j G A^) is a point wise particle, which has its own coordinates f{S), speed v{S), 
and internal shift sh{S). Simultaneously, we will define at each step the new objects 
called bonds. Speaking non formally, a symplex is some set of samples connected 
by bonds with each other, which roll with very hight speed so that their kinetic 
energy is large whereas the speed of the symplex as the whole is low. Factually, the 
swarm will be divided to two fractions: separate samples, travelling with the largest 
possible speed c, and symplexes of orders J > 1, which speed is much lower. This is 
the reason of the itroduction of the time of waiting: this is the time, during which 
the symplex S files, but which is not enough that S can overcome the distance to 
the neigboring cell. As we will see in the next section the fixation of the division of 
space to cells is necessary for the simulation of quantum dynamics; this is why we 
cannot consider the continuous movement of symplexes. 

For this we first of all define the procedure of association of a sample s with 
a symplex S of the order j. Let S be located in the cell f, and s be a sample, 
which has just filed to this cell. We suppose that during the considered step the 
distance between s and Sj reaches its minimal value. We now suppose that the 
time files continuously, it will not change the duration of step, because of the great 
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fraction c/v{S). At the instant when the minimal distance is reached we estabhsh the 
absolutely rigid connection between s and S so these objects start to roll classically 
around their common center of masses, the resulted object we call a symplex of the 
order j + 1. Its speed vi = {v{s)m + v{S)M) / [m + M) where m, M are masses 
of s and S correspondingly, results from the classical consideration of the system 
s + 5", the single non classical object here is the bond we introduced between s and 
S. We call this bond the main bond of the resulted symplex Si and designate it as 
h{Si). We note that due to the absolute rigidity of the bond it does not produce any 
work; the energy, impulse and momentum of impulse are thus conserved in course 
of assotiation. 

The process of bond elimination is exactly reverse of the association; it results 
in the separate sample, which flies away from the symplex that becomes less to one 
in the order. The iteration of the associations q times starting with S gives us the 
symplex S", which order is j + q. Factually, a symplex of the order j represents 
the sequentially nested symplcxcs, each of which is the result of association between 
some sample and the predecessing symplex. The sequential dissociation goes in the 
reverse order. 

Now we are ready to define the step of the main evolution. The step looks as 
the following sequence. 

• For each symplex S we change sh{S) to sh{S)-\-v{S)t and check the inequality 

sh{S) < d where d is the distance to the closest ncigboring cell along the 
direction of vS. If it is satisfied, we go to the next action, if it is violated, 
we shift S to the closest neigboring cell and put its interbal shift to zero: 
sh{S) = 0. 

• Each symplex obtains the additional speed as 6tVV^/m, where V is some ex- 
ternal potential, t is the duration of step, m is the mass of a sample, h is some 
constant. 

• From each symplex S we take [aj] last samples from S and sequentially elim- 
inate bonds connecting them with the symplex, call them belonging to thin 
layer and let them fiy to the closest neigboring cells with the speeds they have, 
where j is the order of this symplex, a is some constant a <S 1. 

• After the arrival of each sample belonging to thin layer to the closest neigboring 
cell we associate it with the corresponding symplex and recalculate its speed 
by the law of association. 

We call the effect of bonds annihilation by peak explosion, because it seems as 
the explosion of the symplex, where the hidden energy of rolling samples transforms 
to the visible kinetic energy of released samples, which fiy out of the cell at hand. 
The third item says that these samples are immediately absorbed by the neigboring 
symplexes; their energy and impulse is thus captivated by these symplexes. Peak 
explosion bears the resemblance with the gas fraction, the movement of symplexes 
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following from the first point we can call the liquid stream. The first and second 
items establish the rules for phase transition: from liquid phase to the gas and vice 
versa. The velocity of gas fraction is much higher than of the liquid phase; this is 
the representation of non relativistic character of qynamics we regard: 

'^symp C. (3) 

It causes the principal difficulty when p becomes too small, because here the corre- 
sponding symplex will obtain the speed comparable with c. In this case our defini- 
tion of symplexes loses its validity at all because this definition presumes that fiying 
samples appears to a cell when a symplex stays; this symples cannot leave the cell 
until the sample associates with it. We note that this difficulty is the same, as in 
Bohm hydrodymanical approach; it is hardly possible to overcome it with swarms. 
On the other hand, the area where p = plays the principal role in the simulation 
of quantum dynamics; its infiuence shows immediately because of a great speed of 
gas fraction. In a quantum swarm, no difference how it is treated (hydrodynamics 
or dynamic diffusion) it is impossible to make consideartions local, as in ordinar 
hydrodynamics when we separate "small cube" and consider a stream through its 
border. The reason is the existence of a gas fraction. Further we will have the more 
quantitative characteristic of this form of quantum non locality. 

We choose the initial state of the swarm in such a form, where each cell is 
occupied by exactly one symplex; the density p(f) will be then proportional to the 
order of the symplex occupying the cell f. 

The evoltion we have defined depends on the constant a, which is the intensity 
of bonds annihilation; we can treat it as the 1/AT where AT is the average time 
of bonds life. The process of bond annihilation will be thus the result of Poisson 
random process with the given intensity. 

We now consider two neigboring cells x and xi with the common border, the 
density of stream through this border Px,xi , and find its variation in the time. Density 
of stream is the stream divided to the square x^ of the border (that is the stream 
through the border of a unit sized cube - we take it to pass from quantity of samples 
to the density), we call it simply the stream. Why we focus at the stream variation, 
not the stream itself? Because the stream depends on the initial state of the swarm, 
which must be given in advance: we cannot derive it from the mechanism of samples 
speedup. The mechanism of dynamic diffusion we defined infiuences to the stream 
only through its variation, which is the dynamical magnitude, as a force. We will 
see, that in the initial moment the stream is created by the movement of symplexes, 
whereas separate samples belonging to thin layer create the variation of stream. 

Let j and ji be the orders of the symplexes in the cells x, and Xi correspondingly. 
We must estimate the deposits to the steam variation from gas and liquid fractions 
separately. 

1). Liquid fraction - inertia. We estimate the deposit of the first item. The 
function of density p typically has the strongly discontinuous form, where "peaks" 
are interspered with "holes" . Of course, it is the infiuence of the fixation of x; further 
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we recognize that it is inavoidable. To cope with this discontinuity we can reform 
the first item, and let samples fly by the annihilation of bonds, as in the second item. 
Here the part of aij samples will fly to the direction / where ai = I ■ v{r), where 
ao is small constant, <^ a. This way of making p more continuous, however, 
contradicts to the invariance towards the change of inertial readout system: our 
mechnism must not depend on this choice. We make this agreement temporarily, 
in order to estimate the deposit of the liquid phase to the change of stream. Thus, 
this deposit will be Vnormir)^ p{r) , where Vnorm is the component of speed of the 
swarm, orthogonal to the border, because it comes from the variation of the order 
of symplexes, which comes through this border. 

2) Liquid fraction - external potential. This deposit equals bpW. 

2) Gas fraction. The deposit of peaks explosion is the In view of the second item 
of our definition of a step, this stream equals the difference between the number of 
samples filed from the cell x to xi, and the number of samples flied in the opposite 
direction, that is the difference between j and ji. The deposit of gas fraction is thus 
aVp. 

We can make the following conclusion. If constants a and b are substantially 
larger than oq, the change of the stream in the unit of time is following equation 

^^ = aVp + bpVV. (4) 

This equation characterizes dynamic diffusion. To show that this mechanism can 
serve as the approximation of quantum unitary dynamics, we have now to derive it 
from Shredinger equation. 



3.2 Reducing of Shredinger equation to swarm 

We now consider the swarm, for which ([2]) is true where is exact solution of 
Shredinger equation, and try to derive the evaluation of the stream variation from 
Shredinger equation directly, accepting some features of sample moving, which agree 
with the mechanism of dynamic diffusion. Our aim is to show that this variation 
satisfies (jlj). 

To ensure the main requirement (I2l) we must define swarm density as 

where is the number of samples occurring in the cube with the side dx with the 
center r, Nfotai is the total number of samples in the swarm. To compare with the 
solution of Shredinger equation we would have to launch in this definition x — > 0, 
which would mean that we consider not one swarm but the sequence of swarms 
with densities pn with increasing n. We will not do it in order to avoid the useless 
complification, instead each time when it is necessary we agree that it is possible to 
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continue our division of the space to cubes so that x will decrease in the admissible 
limits. We write p{x) = |\l'(x)p, which means that 

(n^oo), (6) 

without special mentioning. Such a sequence of swarms realizing the approximation 
to the density of wave function - solution of Shredinger equation we call the admis- 
sible approximation of quantum evolution. Now having this agreement, we fix the 
grain x of spatial resolution. This fixation is necessary for the computing of stream 
variation. We call the stream with the fixed x the detailed stream, to emphasize its 
substantial dependence from x. 

At first we show that there exists the quantum swarm with the local shifts of 
samples, e.g., that the equality ([2]) can be guaranteed by only shifting samples on 
the close distances: between neigboring cells. 

We would easily guarantee the satisfaction of if we do not impose any limita- 
tions on the speeds of samples and on its change. It means that samples can move 
in course of one step on any distance and the equation ([2]) will be true on each step 
with the required accuracy; it will be conserved in some time frame t. However, 
this mechanism is useless because it depends on the knowledge of the wave function 
whereas our aim is to manage without the wave function at all. 

We take up the quantum swarm, and start from Shredinger equation 

d^(r t) h"^ 

^h^^ = -— Avl/(r,t) + V;„i(r,t)v^(r,t), (7) 



which we can rewrite as 



2m J h 



for the real and imaginary parts , \E'* of the wave function = \E''" + z\E'*. We are 
interested in the evolution of density only, e.g., the function 

We fix the value x and apply for the approximation of second derivatives the 
difference scheme of the form 

d'^-^ix) ^'(a; + x) + ^(x-x) -2^(x) . 

(9j 
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for each time instant, provided the wave function satisfies all conditions for such 
approximation. Since the addition of any constant to the potential energy Vpot does 
not influence to the quantum evolution of the density, we can consider instead of 
Vpot the other, equivalent potential Vpot + where a = —j^^q2 that results in the 
dissapearence of the summand 2\E'(a;) in the difference schemes for second derivatives 
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on X, y, z (from that the coefficient 3 appears) after its substitution in Shredinger 
equationj^ For the simphcity of notation we introduce the coefficient 

h 1 

Since we yet do not know the mechanism of moving of samples in the quantum 
swarm, we suppose that we simply take off some quantity of samples from one cell 
or place them to this cell from some storage. We divide the evolution of quantum 
swarm to so small frames of the longitude t, that on each frame samples travel in 
the framework of two neigboring cells. If we prove that the diffusion mechanism 
provides the evolution on each of such frames, it will be true for the whole evolution 
because our supposition about the exchange between two closest cells do not limit 
the generality. We also agree that these cells differ from each other on the shift to 
X along the axis x, which does not limit the generality as well. We denote centers 
of these cells by x and Xi = x + x. 

Now we divide the segment At into two parts: 1) and 2), and agree about the 
following. On the half 1) we will consider only the change of p{x), which is caused 
by p{xi), and on the half 2) - the change of p(xi), caused by p{x) (thhe effect from 
the small change of p{x) on the ffist half has the more high order and we ignore it). 
The density is uniquelly defined by the wave function, hence we will consider the 
change of the wave function on each parts 1) and 2). We show that the required 
change of the wave function can be obtained so: on the first half we consider the 
change of \E' corresponding to the equation 

9^^'(x) ^ ^(x + x) 

^ dx^ ^ (x)2 

and on theh second half - corresponding to the equation 

^^^(xi) ^ ^(xi - x) 
dx"^ (x)2 

Really, if we accept that the wave function changes accordingly to these two equa- 
tions for all nodes, we obtain that the final change resulting from the described 
process on all possible x coinsides with the change resulting from the difference 
scheme (Q. Taking into account that the equation 2) follows from the equation 1) 
after the replacement of x by xi and vice versa, we obtain that the result of the 
sequential actions 1) and 2) can be represented by the system of equations 

vl/[(x) =--f^^{xi) + V{x)^'{x), 

^i(x) = 7^'"(2;i) - 1/(a;)^'"(x), ^ ^ 

^This trick is not accurate from the mathematical view point: we use the fact (possibility to 
add a constant to potential), which is substantiated by analysis, whereas we fix x and in further 
cannot launch it to zero. However, it is not critical: we could preserve the last summand and fulfil 
computations with it; it gives the same result but complicates computations. 
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and the analogous system obtained by the replacement of x by Xi and vice versa, 
where the lower index means the differentiation on t. It means that we regard only 
the exchange of samples through one fixed border between cells x and xi. If we then 
sum the other deposit corresponding to the exchange between x and xq = a; — x, we 
obtain the full change of density within the effect of the separation to two fractions 
(see below). This exchange has the analogous form and all further computations 
can be fulfilled for it with the corresponding result; we thus take up only exchange 
between x and xi. 

The system (fTOj) is true in the supposition that samples move from x to xi. 
When samples move from Xi to x we obtain the analogous system obtained by the 
replacement of x by Xi and vice versa. Shredinger equation will then express the 
result of the general evolution process consisting of the both cases where x and xi 
can settle down by six ways along three coordinate axis. By the time frame At we 
now mean just such short time segment when the exchange goes only between x and 
Xi (and, may be, the storage; we will shortly see that the storage is not needed). 

For such a segment we have 

P.,., = = 2^*(x)(7^'^(xi) - l^(x)^'-(x)) +2^'-(x)(-7^^(xi) - V{x)^'{x)) = 

= 27(^*(x)^'-(xi) - Wix)W^ixi)) = 

(11) 

where by ^^^\x,xi we denote the deposit to the derivative of density on the time, 
which appears from the moving of samples through the border separating cubes that 
contain points x and Xi. It gives that the decrease of samples in one cell equals the 
increase of them in the other, e.g., the evolution of quantum swarm satisfies the 
condition of locality and we can speak about the stream of samples through the 
element of surface when we fix the value of x. 

3.3 Dependence of swarm dynamics from grain 

We approach to the most important thing in the description of quantum dynamics: 
the dependence on the grain x. We cannot confirm that the expression ffTTj) deter- 
mines the number of samples, which increases the density of the corresponding cell 
when they pass through the considered surface in the unit of time. Expression for a 
detailed stream (fTTj) is not expression for an ordinary stream of particles through a 
surface. The reason in that it critically depends on a value x, and it is impossible to 
direct this value to zero. We cannot apply to a swarm the fundamental reception of 
the mathematical analysis, consisting that it is possible to divide unlimitly intervals 
{x,xi) and to pass to a limit at x — > so all values of group magnitudes (the 
density and speed) will keep the values (and even them will specify). Expression 
f lTTj) indirectly specifies that at the swarm there are samples with essentially differ- 
ent speeds so we cannot speak about speed of all swarm in the given point. Let's 
recollect that the swarm has two parts: fast (separate samples from the thin layer) 
and slow (symplexes). The mechanism of moving of samples of the fast part consists 
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that each of them on each step At jumps in a direction of its speed on the number 
of cells proportional At. The mechanism of moving of a sample of the slow part of 
the swarm consists that it waits number of the steps, proportional 1/At, being at a 
stop, and only then moves exactly on one step in a direction of the speed. In this 
case we cannot write for the detailed vector stream (density) p expression 



dp{r,t) 1 



dt A^(x)3 

S{t) 



p{r,t)n{ri)ds{ri) (12) 



that is true for the ordinar stream p. Instead we must use the expression 

^^^"■'^^ - ^ I' p{r,t)\(^,^,,jn{n)ds{n) + A, (13) 



dt N{xy 

S{r) 

where the summand A have the following nature. 

Let X and At be fixed. To write expression ( fT2l) . it is necessary to guarantee 
that all samples of the fast part of swarm, which have past through the border in 
a direction r, remain in the corresponding cell, and do not jump out of it because 
of the great speed. For this, it is necessary to make At small enough. But then 
samples of the slow part of the swarm will not have time to move a little at all! 
For maintenance of movement both fast, and slow samples there is only one way: 
to reduce a grain x. But we also cannot make it, because then at once will increase 
speeds of the fast part of the swarm, and we should consider already the other 
swarm. 

Thus, it is impossible to replace discrete character of moving of samples of fast 
and slow parts of the swarm with continuous movement, as in case of the classical 
continuous environment. Summand A it makes sense the phase transition between 
fast and slow parts of the swarm. We turn at the necessity to consider two fractions: 
slow liquid and fast gas. (Otherwise, if we make possible for symplexes move at each 
step, we would have to consider jumps of a thin layer samples: the jump of a fast 
sample means that it disappears in one cell then appears in other cell which can 
have no common borders with the first.) There is no the usual stream in a quantum 
swarm, we can thus operate with the detailed stream only. 

Inapplicability of usual analytical receptions to the quantum swarm does not 
mean impossibility of a substantiation of a way of modelling offered by us by means 
of classical mathematics at alll I do not exclude a finding of such substantiation, 
but - only for the case of one real particle. However, it is impossible to find a classical 
substantiation of the method of collective behaviour in its general form, for many 
particles, because here the limitation of quantity of samples plays the central role 
that makes impossible application of the ideology of the mathematical analysis. It, 



^In a case the delta of function of Dirac the classical mathematical substantiation formally has 
been found in the form of linear functionals. 
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in particular, means also hopelessness of attempts to prove formally that our mecha- 
nism of moving of samples of the swarm in all cases gives good approximation of the 
exact solution of Shredinger equation. Therefore the following below substantiation 
of the proposed mechanism should be treated as the explanatory to algorithm work, 
but not as the strict proof of its suitability in all cases. The general ideology of 
constructivism (see [6j) offers us one: to rely on direct computer modelling. 



3.4 Reduction of Shredinger equation to dynamic diffusion 

We introduce the following parameters depending on x: 

/ = — ^ K = -^ (14) 

where h is Plank constant. The parameter / is called the intensity of action of the 
density gradient, k is the intensity of potential. 

Now for construction of the mechanism of change of speed of samples we will 
find the change of the detailed stream ■^p\{x,xi) ndS of quantum swarm through a 
surface of a small cube, the normal to which is parallel to the axis OX. For this 
purpose it is necessary to differentiate the expression (fTTj) on the time: 

p'^ = 27[(7^"(xi) - V{x)W'{x))^''{xi) + ^*(x)(-7^*(x) + V{xi)^'{xi))- 
(-7^^(xi) + 1/(x)^*(x))^*(xi) - ^'■(x)(7^'^(x) - V{xi)'^^{xi))] = 
272(^^(xi))2 - 27l/(x)^"(x)^"(xi) - 2-f\^'{x)y+ 
2jV{xi)'^'lx)^\xi) + 2^^{^'{xi)y - 27r(x)^^(x)^^(xi)- 
272(^^(x))2 + 27\/(a;i)^''(x)^^(xi) = 
272((^-(xi))2 + (^*(xi))2 - ((^^(x))2 + (^H^))')) + 
2'y[{V{xi) - V{xM^''{x)f + {^\x)f) + o(x)], 

(15) 

where o(x) = (^^(x)^^(xi) + ^^(x)^*(xi) - ((^'■(x))^ + {^'{x)y)){V{xi) - V{x)). 
We thus can write for the detailed stream the formula 

Ap| = {-IV p - KpVV)At. (16) 

We see that Shredinger equation gives the same change of stream as the dynamic 
diffusion (see equation (jl])) if we put a = —I, b = —k. 



3.5 Restoration of wave function from dynamical diffusion 
swarm 

For a conformity finding between the standard description of a quantum state 
through its wave function \1/ and its swarm representation, we will assume that 
the last is described by pair of functions 

p(t,f), p{t,f), (17) 
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where the scalar function p is the density of samples, and the vector of impulse of 
the swarm p{r) is the sum of speeds of samples occurred in the small cell with the 
side X and the center in the point r. The swarm impulse does not include weight of 
the real particle; p(r) there is the swarm characteristic in which "the weight" role 
plays the total number of samples. 

Such pair does not use concept of complex number, and does not give the beau- 
tiful differential equations of Shredinger type for p and p. Moreover, the mechanism 
of dynamic diffusion introduced by us for imitation of quantum evolution, consid- 
erably differs from classical processes (for example, a heat transfer or fluctuations) 
that its intensity depends on the chosen grain of the spatial resolution. We have 
gone on it for the sake of the main thing: economy of computing resources which 
are not simply necessary for modelling of quantum dynamics of complex systems, 
but are absolutely necessary for its research in general. 

Given a state of swarm we will show, how the usual complex wave function ^! 
can be restored from the pair f|T7|) . We accept that the carrier of wave function is 
connected, that is any two points in it can be connected by a curve which is not 
crossing area of zero density p = 0. This restriction is related to that is available 
in the diffusion Monte Karlo method. For this purpose we consider the equality 
( ITTj) . and substitute in it the expression for wave function through density: \E'(r) = 
a/ p(r) exp(z0(r)). The problem consists in the calculation of the phase 0(r) of 
wave function. We will notice that as the relative phase between various points has 
physical sense only, we can fix some point r and consider a phase of other point ri 
relatively to r. If ri is close to r, the equation ffTTj) gives us 



(r) — 0(ri) = arcsin k{x) 



2 P(r-ri) 



^/p{r)p{ri) 

that leads to the following formula for the relative phase: 

0(^1) = / k{x)^v d'y 



where the path 7 goes from r to ri. This equation explicitly depends on the choice 
of the path 7 hence we have to prove its correctness, e.g., independency from the 
choice of 7. 

We mention that this derivation will be correct only in the case when p > e > 
for some constant e > 0, e.g., the density must be separated from zero in all area of 
consideration. Since the phase is determined only within to the integer multiple of 
27r, different choices of the path can lead at most to the choice of the phase to such 
a number that takes place, for example, for excited states of an electron in hydrogen 
atom with nonzero magnetic number. We show that the integration of the speed v 
of swarm along a closed path preserves its value in the time the more exactly, the 
lesser grain of spatial resolution x is. It results that if in the initial time instant 
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the definition (fTSi) is correct it preserves the correctness for the following instants 
as well. 

We thus consider the derivative of the integral of the speed of swarm along the 
closed path 7c. Applying the formula ( 1T6|) and taking into account that dp/dt is 
proportional to p dv/dt, we obtain 



for some A, B. The first summand gives zero after the integration along the closed 
path because it is V In p, the second summand gives zero by the analogous reason. 

It is now sufficient to check that the definition f lTSj) is correct in the initial instant 
that can be done straightforwardly for each task. In the case when the wave function 
of initial state can be obtained from the ground state of electron in hydrogen atom, 
where w = 0, the correctness follows from the proven because there is no phase shift 
to 2'Kk. The phase of each groundstate does not depend on the point. If for the 
obtaining of the initial state in the considered problem we have to start from some 
escited state we must at first check the correctness for this state. 

On the basement of these computations we can write formulas connecting swarm 
parameters with the wave function: 



for some a, h. These formulas permits to pass from the wave function to the swarm 
and vice versa. This description has two features. At first, beyond Shredinger dy- 
namics on the level of grain x stays the swarm dynamics of the lower level with grains 
Ax, At, so that there is the substantial dependence of swarm parameters (intensity) 
from the chosen grain of spatial resolution x. At second, quantum dynamics pre- 
supposed the presence of the thin layer and the specific behavior of samples towards 
it. 

In formulas f l2U]) the mean speed v of swarm approximately equals the mean 
speed of such a part of swarm, which does not contain thin layer in each point 
where the density does not converge to zero. There the thin layer does not give 
the substantial deposit to the impulse of swarm because in each small cube there is 
only one sample from thin layer. If we apply to the quantum particle the external 
potential Vpoti it brings the deposit to the impulse of swarm because it infiuence to 
the main part of swarm that does not belong to the thin layer. 

4 Case of many particles 

Now we show that the method of dynamic diffusion can be generalized of the case of 
many quantum particles. Let we be given the set of n quantum particles, which we 




(19) 




(20) 



V = bV4>{r), 
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enumerate by natural numbers: 1, 2, . . . , n. To find the efficient scheme of modelling 
algorithm we should apply sequentially the method of collective behaviour, at which 
the algorithmic reduction of quantum states is the inbuilt property. The most correct 
decoherence model, which we named absolute, says that decoherence is the reduction 
of quantum state 

i^) = E^^b') (21) 

j 

while the memory of the modelling computer cannot contain full record of this 
state. As we have shown above, such model gives Born's rule for calculation of 
probabilities by an outcome of measurement of a quantum state of system, and 
it proves a correctness of the given model. But such form of absolute model of 
decoherence cannot yet serve as heuristics for modelling algorithm as we while do 
not have any way of modelling of unitary quantum dynamics for many particles, 
except calculations within the limits of matrix algebra, and this way as we saw, is 
too expensive. 

We will consider swarm representations of ours n particles 1,2, ...,n where 
Si, S2, . ■ . , Sn are their swarms, corresponding to their states 

1^1), 1^1/2), • • • , l^n)- If to consider the ensemble consisting of all these samples, it 
will represent a simple state of a kind 

1^1) (S^ 1^2) (S) • • • l^n)- -B^t for representation of the entangled state of a kind 

*)= Yl ^ilJ2,...jJjl,J2,...,Jn) (22) 

jl,j2,--;jn 

we need to introduce a new essential element into the method of collective behaviour. 
These are bonds between samples of different swarms, which we call own bonds, to 
distinguish them from non own bonds, which we have introduced earlier to join sam- 
ples in symplexes. The basic state ji can be considered as coordinates of the particle 
i in the corresponding configuration space. Representation of the wave function in 
the form (!22|) means that there are bonds, which connect points ji, j2, • • • , jn in one 
cortege. 

In the method of collective behaviour (see [^) we will accept that bonds connect 
not spatial points, but samples of various real particles. These bonds can be written 
down as corteges 

S = {Si,S2,...,Sn) (23) 

where for any j = 1,2, ...,n Sj G Sj. Wave function |$) then is represented as 
set 5* of corteges s so that for any j = 1,2, ... ,n and Sj G 5*-,- there is exactly one 
cortege of the form fl23|) . Each cortege plays a role of the so-called world at many 
world interpretation of quantum theory. We consider this cortege fl23l) as a sample 
of system of n particles. In case of one particle we saw that action on a thin layer of 
the force proportional to the gradient of density of the swarm can simulate quantum 
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dynamics. We will see that this process is directly generalised also on the case n 
particles. We name S swarm representation of the system of n particles. 
The swarm density S is defined as 

^Jio;^^%^, (24, 

where iVri,r2,. dx is a total number of the corteges which have appeared in the 3n 
dimensional cube with the side dx and the centre ri, r2, . . . , r„, N - total number of 
corteges. 

If wave function |$) is the tensor product of one partial wave functions: 



i=l 

the corresponding bonds then can be obtained by the random choice of samples 
from the uniform distribution Sj G 5*^ for everyone j = 1,2, ... ,n which will thus 
form each cortege Si, S2, ■ ■ ■ , s„. With such choice of corteges we will receive that 
the density of a corresponding swarm satisfies to Born's condition, which can for 
swarms be written down as 



m\'-^ (25) 



where D C i?^", N^^§ is the total number of corteges which have appeared in area 
D. But for the entangled state |$) such choice of corteges for a set of swarms S 
will not give us a condition fl2S]) . We thus should take ( HM for definition of a choice 
of corteges in S. But for definition of a swarm we should define also speeds of all 
samples, namely, generalise equahty (!20|) on the case of n real quantum particles. 

Let \l/(ri, r2, . . . , r^) be a wave function of system of n particles, 
\1/ = |\l/|exj9(i0(ri, r2, . . . , r„)) its Euler decomposition. We will designate through 
Vj0(ri, r2, . . . , r„) the gradient taken on coordinates of the particle j where 
j G {1,2, ... ,n} there is a fixed number. Generalisation of formulas (120 p on n 
particles looks like 

|^(f)| = 

0(r) = / kv d-f, ('26) 

7: ro — >r 

V = aV0(f), 

where f designates ri, r2, . . . , r„, V designates (Vi, V2, . . . , V„), and 7 is a path in 
3n dimensional space. ( l26l) is enough rule for the definition of the swarm on the 
given wave function if we agree to unite samples in corteges irrespective of their 
speeds. At transition from the case of one particle to the case of many particles it 
is necessary instead of a sample of the swarm for one particle everywhere to insert 
the sample of the swarm for all system of many particles. 
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The external description of change of speeds will be the direct generalization of 
the case of one particle. The gradient Vp of the density of this swarm will cause 
change of speed of each sample of a separate particle in those corteges, which belong 
to the thin layer. Return of an impulse of the thin layer to all swarm occurs in 
each cell of configuration space in the same way, as in case of one particle. Dynamic 
smoothing of formed peaks also goes by the general rule. 

Let's consider now the swarm of samples of quantum system with n particles. 
Each of them represents Everett quantum world, and looks like s = (si, S2, ■ ■ ■ , Sn)- 
How such cortege receives an increment of speeds Aps = {Api, Ap2, . . . , Ap„)? For 
this purpose it should be in the thin layer. Let it already there is. Then we will 
consider all other corteges lying close to s in sense of the metrics of n particle 
configuration space, and we will find the gradient of their density which is looking 
like: 



Now each sample Sj receives an increment of its own impulse on —IVj. Cortege 
elements (|27|) are not defined by density of separate particles. They are defined 
by density of the swarm for all system with n particles. That is influence on Vj 
renders not only swarm density of j-th particle taken separately, but samples of 
all other particles which are connected in one cortege with an environment of j-th 
particle provided that these corteges lie close to s. We could consider only density 
of separate particles only in the case when the state is not entang ledl Thus if we 
wish to study the general case of the entangled system in swarm representation, 
we should enter into consideration bonds between samples in individual swarms, in 
particular, for corteges belonging to a thin layer. As in the case of one particle, we 
redefine thin layer at each step. 

Here again in case of a non entangled state the role will play only affinity of 
samples of separate particles. Bonds between the different samples entering into 
one cortege will be involved only if we have an entangled state. Thus, the thin layer 
in a case of the entangled state of n particles also cannot be received as a simple 
combination something like "thin layers" for separate particles. It is defined through 
bonds in corteges, that is has essentially multipartial character. 

Let m be the number of samples of each separate particle. Then, owing to our 
condition of uniqueness of the cortege containing the given sample of any particle, 
we have exactly m various corteges, which are not overlapping pair wise If we launch 
the number of real particles n to infinity, having left m to constant, we with growth 
of number of real particles n will receive the increasing deviation of our model 
from the exact solution of the corresponding Shredinger equation. Thus, fixing m 
means presence of the internal factor of decoherence, which is independent from 
any "environment". On the other hand, complexity of model will linearly increase 
with growth n, instead of exponential, as in the standard formalism. That is the 
method of collective behaviour realises absolute model of decoherence which in it 

■^just for non entangled states the classical consideration of ensembles is legal. 




(27) 
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comes owing to limitation of memory of the modelling computer. 

Acceleration of a thin layer as a result of action of a gradient of density of 
a swarm in case of limitation of number of samples gets the specific form. If to 
divide configuration space of each real particle on s cubes the configuration space 
for n particles will appear divided on s"" cubes. Therefore we can speak about any 
"calculation of a gradient of density" in Bohm' sense only in that case, when m ^ 
that is when it is not a lot of real particles. If this inequality is not true, division of 
the force operating on a separate cortege from the swarm on " analytically described 
component" and "processing of peaks" loses meaning. "Peaks" will be actually 
everywhere, and it is necessary to describe process of their dynamic smoothing 
correctly. 

Now we take up the mechanism of changing of speeds for samples of n particle 
system. It is formulated as in the case of one particle, but the close samples will 
be replaced by close corteges. Namely, a cortege we treat as a symplex of the first 
order. If two corteges Si and ^2 are close in the sense of n particle configuration 
space, e.g. they belong to the same cell in this space. Then in the instant when the 
distance between then is minimal the own bonds between the corresponding samples 
are established and we obtain the n particle symplex of the second order. We call 
the sequence of these one particle bonds an n particle bond, where one particle 
bonds arc its components. A symplex for the n-th order arises from a symplex of 
n — 1-th order and a sample by the same procedure. The decay of all one particle 
components of a bond happens simultaneously, that is we can speak about the decay 
of n particle bond in a given time instant. The decay results in the flying away a 
sample, etc. Thin layer consists of samples, which move from cell to cell at the 
short time frame. With this mechanism we again obtain the action of gradient of n 
particle swarm of corteges, which acts on the thin layer. 

5 Computer simulation of dynamic diffusion 

This computer model should be considered as a proof-of-concept of Dynamic Diffu- 
sion Swarm. 

Thus, was decided to model few base experiments inside 1-dimensional rectan- 
gular potential hole: dispersion of an gaussian wave-packet, ground state of wave- 
function (in comparison with linear modulus function), collision of two gaussian 
wave-packets (as illustration of quantum interference effects) and oscillation of a 
particle between two potential holes, separated by high potential barrier. 

Algorithm's convergence is very close to Implicit Runge-Kutta after a little time 
scheme collapses, but consistently shows quantum effects. 
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5.1 Algorithm overview 



Workflow of algorithm is an evolution iterations of simplexes array and it's inter- 
action with potential barriers array (potential hole is represented as two very high 
potential barriers). Scheme of simplex creation will be provided in this way: 

SimplexName ( speed, samples, position) 

. There are three phases: Explosion, Flight, Rearrangement. 

• Explosion. Coefficient a indicates amount of samples ([0; 1]) detached from 
a simplex (see Mechanism of bond creation and annihilation). In place of 
simplex, three others are created: one will then flight to the left, other one 
to the right and third remains in place. This is a pseudocode of Explosion 
phase (let original simplex be Original(original.samples, original.speed, origi- 
nal.position)): 

for each original simplex in array 

create new simplex with speed: (original.speed - max_speed) 
samples: (a * original . samples ) 
position: (original . position) 
create new simplex with speed: (original.speed + max_speed) 
samples: (a * original . samples ) 
position: (original . position) 
create new simplex with speed : (original . speed) 
samples :( (l-2*a) * original . samples ) 
position: (original .position) 
remove original simplex from array 

• Flight. Here is a pseudocode of phase: 

foreach simplex and foreach barrier 
if simplex flights through barrier 
handle barrier 

simplex . position = simplex . position + simplex . speed 

• Rearrangement. As simplexes flight continues, and do not interact with each 
other, rearrangement phase is needed to perform interaction of nearby sim- 
plexes and perform some kind of measurement, arranging simplexes back to 
grain-dependent grid. Here is pseudocode of rearrangement phase: 

foreach point on grid 
find all nearby simplexes 
interact simplexes 
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5.2 Results 



Comparing ground state with artificial linear modulus function. Systems 
evolve for the same time. On screenshots it is obvious, that ground state 
appears to be more stable: 




(a) t = (b) t = 1800 

Figure 1: Ground state 




(a) t = (b) t = 1800 

Figure 2: Modulus function 



Numeric comparison of convergence was analyzed with a metric: 



TV 



M {wave function, T) = 'y~^^\wavef unction. samplesl.J^— wave function. samples 



t=T\ 
x=i I 



1=1 



Result is present on next figure: 



- M(niadjlus} — lvl(ground state) 



Figure 3: M{modulus, t) - M(ground state, t) 
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• Interference of two wave-packets. Results were compared to Schroedinger's 
equation numeric solution, performed in Wolfram Mathematica. Firstly, this 
experiment was modeled in Wolfram Mathematica with implicit R-K, explicit 
R-K and 'Automatic' computation method. Implicit method is absolutely in- 
consistent. Explicit method does not show interference effect and disconverges. 
'Automatic' results are represented on these two figures: 
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Figure 4: Two gauss wave-packets, Wolfram Mathematica 
And this are results, received in DDS modeling: 




(a) t = (b) t = 2100 

Figure 5: Two gauss wave-packets, DDS 



• Oscillating of a particle between two potential holes, separated by high po- 
tential barrier. Firstly, should be marked, that ground state disperses faster, 
than cycle of oscillating occurred. Thus, was applied an artificial smoothing 
transformations. Stability of ground state, represented by DDS is a theme for 
further investigations. Results with smooth-out transformation: 

5.3 Modeling conclusion 

These results can be interpreted as a proof of constructive approach to wave-function 
dynamic modeling. As was described above, this model can be expanded to multi- 
particle case, which is impossible with classic Schredinger's equation. Certainly there 
is a lot of work to be done, before it can be successfully used in real computations 
of many particles systems. 
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Program built on .Net WPF platform. Model is architectured the way sepa- 
rating Model and View, so it can be integrated into other software (for instance, 
Experiment Work) . 

6 Conclusion 

We show how quantum unitary dynamics can be simulated in terms of a swarm of 
classical point wise samples, that is the embodyment of prequantum dynamics in the 
sense of [5]. Mechanism of speedup samples rests on the new object - bonds, which 
are absolutely rigid connections joining samples into symplexes. A symplex is a set of 
close samples, which roll with high speed, whereas a symplex as a whole moves with 
low speed. Quantum dynamics comes from the creation and annihilation of bonds 
that formed a kind of Poisson process. In course of it samples fly out from their 
parent symplexes and rapidly jump to the other close symplexes, giving them their 
impulses. This process results in the rapid transmission of impulse. The derivative 
of a stream of samples through the small border between cells turns proportional to 
the density gradient, where the coefficient depends on the grain of spatial resolution 
as x~^. 

This dependence makes impossible to launch the grain to zero and hence, to use 
the exact analytic apparatus to the swarm dynamics. The reason of it lies in the 
nature of diffusion swarm - it consists of two fractions: slow (symplexes) and fast 
(separate samples, fling with the speed of light). The adventage of dynamic diffusion 
is in that this approach to simulation of quantum dynamics can be easily generalized 
to the many particle case, because the mechanism of samples speedup depends only 
on the density in the current point, but not on the closest vicinity of it. In this 
generalization a sample of the whole n particle system is the cortege, consisting of 
samples of the separate particles. This is the idea of collective behavior for swarm 
quantum dynamics. These corteges do not intersect and their total number equals 
the total number of samples of one particle. The limitation of this number plays the 
role of decoherence when the number of real particles converges to infinity. In this 
case the approximation via collective behavior becomes farer and farer from unitary 
quantum dynamics. 

Rotation of samples inside of symplexes bears a resemblance to the spin of par- 
ticle (see [1]). We cannot claim that there is the conventional quantum mechanical 
spin, this spin we can call "hydrodynamic spin". The establishing of the relation of 
this "spin" with the usual spin represents the interesting problem. 

The check of validity of dynamic diffusion is the subject of computer simulation. 
The main problem arises with the area where p ~ 0. For example, for the particle in 
the rectangular potential hole its ground state has the form sin^ax. The influence 
if a wall we represent as the elastic reflection of samples and symplexes from the 
walls. The density then will have peaks along all area, and they move the faster the 
less is the density. Just the existence of these peaks preserves the form of siri^ for 
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the sufficiently long time in comparison to the other functions. The other example 
is oscillation of a particle between two potential holes separated by high potential 
barrier. The huge velosity of samples in thin layer makes possible to overcome this 
barrier, if the general speed created by slow symplexes has the same direction. If 
the direction of this slow part of a swarm is reverse, thin layer cannot overcome the 
barrier. This effect gives the known oscillation, following from Shredinger equation. 

The further development of collective behavior lies to the effect of decoherence 
in complex systems (see [9], [8]). Here the approximate representation of quan- 
tum dynamics in the form of dynamic diffusion will have the adventage over Bohm 
approach, because it does not require the smooth shape of the density. 
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